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ao) a sat ae Background: In South Korea, an outbreak of Middle East respiratory syndrome (MERS) 

University, 77 Jeongneung-to, occurred in 2015. It was the second largest MERS outbreak. As a result of the outbreak 

Seongbuk-gu, 136-702 Seoul, in South Korea, 186 infections were reported, and 36 patients died. At least 16,693 peo- 

pee ple were isolated with suspicious symptorns. This paper estimates the basic reproduc- 
tion number of the MERS coronavirus (CoV), using data on the spread of MERS in South 
Korea. 


Methods: The basic reproduction number of an epidemic is defined as the average 
number of secondary cases that an infected subject produces over its infectious period 
in a susceptible and uninfected population. To estimate the basic reproduction num- 
ber of the MERS-CoV, we employ data from the 2015 South Korea MERS outbreak and 
the susceptible-infected-removed (SIR) model, a mathematical model that uses a set of 
ordinary differential equations (ODEs). 


Results: We fit the model to the epidemic data of the South Korea outbreak minimiz- 
ing the sum of the squared errors to identify model parameters. Also we derive the 
basic reproductive number as the terms of the parameters of the SIR model. Then we 
determine the basic reproduction number of the MERS-CoV in South Korea in 2015 as 
8.0977. It is worth comparing with the basic reproductive number of the 2014 Ebola 
outbreak in West Africa including Guinea, Sierra Leone, and Liberia, which had values of 
1.5-2.5. 


Conclusions: There was no intervention to control the infection in the early phase of 
the outbreak, thus the data used here provide the best conditions to evaluate the epi- 
demic characteristics of MERS, such as the basic reproduction number. An evaluation 
of basic reproduction number using epidemic data could be problematic if there are 
stochastic fluctuations in the early phase of the outbreak, or if the report is not accu- 
rate and there is bias in the data. Such problems are not relevant to this study because 
the data used here were precisely reported and verified by Korea Hospital Association. 


Keywords: Basic reproduction number, Susceptible-infected-removed (SIR) model, 
Middle East respiratory syndrome (MERS), MERS coronavirus (MERS-CoV) 
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Background 

The Middle East respiratory syndrome (MERS) is caused by a coronavirus (CoV), the 
MERS-CovV. In Saudi Arabia, the first case of the disease was reported in 2012 [1]. The 
first case of MERS in the Republic of Korea was identified on 20 May 2015 [2]. A sig- 
nificant outbreak of MERS occurred in South Korea and lasted for almost three months, 
from May to July 2015. 

The 2015 MERS spread in South Korea is the second largest outbreak recorded to date 
[3]. As a result of the outbreak in South Korea, 186 infections were reported, and 36 
patients died. At least 16,693 people were isolated with suspicious symptoms [4]. This 
paper evaluates the basic reproduction number of MERS-CoV, using data from the 2015 
South Korea outbreak. 

The basic reproduction number (generally denoted as Ro) of an epidemic is defined 
as the average number of secondary cases that an infected subject produces over its 
infectious period in a susceptible and uninfected population [5-11]. It can estimate the 
growth rate of an infectious disease at the early stage of the outbreak, when most indi- 
viduals are susceptible [12]. The basic reproduction number of an epidemic is useful for 
determining whether an outbreak of the disease will occur or not [13], and for analyzing 
epidemic properties of the disease further [14]. 

Note that Ro is referred to as the basic reproductive (or reproduction) number (or 
ratio). The basic reproductive (or reproduction) rate is incorrect nomenclature because 
Ro is a dimensionless number that is not related to any physical quantity corresponding 
to rate. 

Based on [11], the reasons to estimate the basic reproduction number of an epidemic 
are summarized as follows: First, we can relatively evaluate the risk of the correspond- 
ing epidemic using Ro. In other words, we can compare the infectivity of the epidemic 
with others, already familiar to us. Second, the reproduction number can be evaluated 
multiple times, e.g., before (and after) an infection control measure intervention. To this 
end, it is needed to distinguish the reproduction number after control intervention from 
the basic reproductive number Ro, which is estimated before the intervention. Now we 
refer to the reproduction number after control intervention as the effective reproduction 
number, denoted by Reg. Then we can compare Reg with Ro, and we can evaluate the 
efficacy of a control measure quantitatively based on Reg. By doing so, eventually we can 
determine how to apply control intervention to reduce Reg to less than one. If Reg < 1, 
it is concluded that the control intervention works effectively, and that the outbreak will 
eventually be controlled by reducing the reproduction number to less than the threshold 
level, 1. 

In this paper, to estimate the basic reproduction number of the MERS-CoV, we employ 
data from the 2015 South Korea MERS outbreak and the susceptible-infected-removed 
(SIR) model [13, 15, 16], a mathematical model that uses a set of ordinary differential 
equations (ODEs). Because the availability of epidemic data is limited, we usually employ 
non-structured deterministic models to evaluate Ro [11]. 

Based on the data reported from the 2015 South Korea outbreak of MERS, we evalu- 
ate the basic reproductive number of the virus, MERS-CoV. We fit the model to epi- 
demic data from the South Korea outbreak, and identify model parameters and the basic 
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reproduction number, Ro. Note that other epidemiological parameters, such as incuba- 
tion period and serial interval, have been discussed in [17] for the outbreak. 

Preliminary work relating to this paper is presented in [18]. This paper includes an 
analysis of the derivation of the basic reproduction number, careful screening of the 
reported data, and a sophisticated approach using the sum of the squared errors to eval- 
uate the basic reproduction number precisely. 

A number of papers and books have been dedicated to the study of Ro for other infec- 
tious diseases. A few of them are as follows: see [19] for the Ro of severe acute respira- 
tory syndrome(SARS), [20] for the Ro of influenza, [21, 22] for the Ro of Ebola, and [23] 
for the Ro of malaria. 

Although some literature on the study of Ro of the MERS-CoV has been reported 
such as [12, 24], the MERS outbreak in South Korea is unique [25]. MERS spread almost 
naturally without any intervention in the early stage, and the Korean government did 
not respond appropriately [3, 25]. The list of medical facilities involved was not even 
announced to public. Ironically, it is the ideal condition to fit a mathematical model to 
the clinical epidemic data and to evaluate epidemic properties of the MERS-CoV, includ- 
ing the basic reproduction number. 

The paper is organized as follows: In “Methods” section, a mathematical model, 
which comprises a set of ordinary differential equations, is introduced, and an estima- 
tion method is discussed for the basic reproductive number. “Results” section describes 
the evaluation result, providing further discussion. “Conclusions” section concludes the 
paper by suggesting additional research. 


Methods 
In this section, we briefly discuss the method employed in this paper, including a defini- 
tion of the basic reproduction number, and we introduce the SIR model. 


Basic reproduction number 
The basic reproduction number is defined as the number of secondary cases that one 
infected primary subject causes on average in an uninfected and totally susceptible pop- 
ulation, over the infectious period [8—10]. Based on this definition, we obtain a math- 
ematical description of Ro via the so-called ‘survival function’ [7, 11]. 

Considering a large population, this description is given by 


Ro = | b(a)F (a)da (1) 
0 


where F(a) is the survival function describing the probability of a newly infected patient 
will be infectious at least to time a, and b(a) is the average number of infected subjects 
whom one infected patient will produce per unit time at time a [7, 11]. The notation of 
(1) follows the usage therein. 

Formula (1) is derived from the Ro definition and can thus be used for any mathemati- 
cal model, not just models given by ODEs. However, it requires explicit expressions for 
F(-) and b(-), which are functions of time. This paper employs the SIR model described 
by ODEs that is introduced in the next section. Because the model does not provide 
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explicit descriptions for F(-) and b(-), we use an alternative expression for Ro, derived 
from the SIR model. 


SIR model 
To investigate the spread dynamics of epidemics, several nonlinear mathematical models 
have been studied. We employ one of the models, the SIR model that is described by 


S$ = —ASI, 
I= BSI —vI, (2) 
R= vi, 


where the states S, J, and R correspond to the number of susceptible people, the number 
of infected, and the number of removed, respectively. Note that state R includes both 
deceased and recovered patients. For each subject group of S, J, and R, we assume that 
the properties of the subjects (for example, infectiveness, susceptibility, and so forth) are 
homogeneous. The parameter is the disease transmission rate and the parameter v is 
the removed rate. Note that both parameters are positive real. 

The notation of (2) follows that used in [16]. Table 1 presents descriptions of system 
parameters and state variables for model (2). Parameters « and t are discussed further in 
the “Results” section. For a further explanation of the SIR model, see [13, 15, 16]. 

Any term related to birth and death in the population that is not caused by MERS 
is not included in the model (2). The dynamics of the disease (e.g., infection or recov- 
ery) is assumed to be significantly faster than that of birth and death in the population. 
Generally, epidemic models such as SIR do not include birth and death because zero 
net change of the population is assumed. If we model an infectious disease with com- 
paratively slow dynamics (e.g., an endemic disease), we must consider dynamic terms 
describing birth and death. 

The system parameters, which are all rates, are positive and real in the model (2). 
Because of the definition of Ro, an alternative description for Ro can be derived from 
model (2): 


Ro = xtd, (3) 


Table 1 Descriptions of system parameters and state variables of model (2) and the equa- 


tion (3) 

State variables, parameters Descriptions 

S umber of susceptible subjects 

| umber of infected subjects 

R umber of removed subjects 

B Disease transmission rate 

v Removed rate 

T Transmissibility of the infection 

K umber of transmittable contacts 


by infected patient per unit 
time 
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where d is the infection duration, « is the transmittable contact number by one infected 
subject per one unit time, and T is the transmissibility of the infectious disease, which 
corresponds to the probability of infection per one contact between an infected patient 
and a susceptible individual. 

It is notable that the removed rate is reciprocal to the infection duration by the 


assumption of constant rates for the SIR model (2): 


oe 
v 
Results 
In Table 2, the history of the MERS-CoV spread status is presented. The Ministry of 
Health and Welfare, Korea officially announced the data. “Infected” represents the 
accumulated number of infected patients. “Deceased” is the number of dead subjects. 
“Recovered” is the number of individuals returning to healthy status. All entries in the 
table are as of the “Date” The number of infected patients includes both removed and 
recovered patients. 

Table 2 shows the MERS-CoV spread data for only the initial phase of the outbreak, 
ie., from 20 May 2015 to 12 June 2015. 

On 7 June 2015, the South Korean government disclosed to the public the list of all 
hospitals exposed to MERS-CoV, with the dates and duration of exposure [4]. This is 


Table 2 Accumulated MERS-CoV patients in Korea, 2015 [4] 


Date Infected Deceased Recovered 
20 May 2 0 0 
21 May 3 0 0 
22 May 3 0 0 
23 May 3 0 0 
24 May 3 0 0 
25 May 3 0 0 
26 May 5 0 0 
27 May 5 0 0 
28 May 7 0 0 
29 May 13 0 0 
30 May 15 0 0 
31 May 18 0 0 
| June 25 1 0 
2 June 30 1 0 
3 June 30 3 0 
4 June 36 4 0 
5 June 42 > ] 
6 June 64 5 ] 
7 June 87 5 ] 
8 June 95 7 2 
9 June 108 7 3 
10 June 122 9 4 
11 June 126 10 7 
12 June 138 13: 9 


Chang BioMed Eng OnLine (2017) 16:79 Page 6 of 11 


the first intervention of the government to control the spread. Before this date, there 
was no control action that could affect estimation of the basic reproduction number of 
MERS-CoV. 

The incubation period of MERS-CoV that can range from 2 to 14 days, is 5 days on 
average [26]. Thus, we use reported data up to 12 June 2015. 

The total population size is denoted by N(Z): 


N(t) := S(t) + I(t) + R(f). 
It can be seen that 


N(t) = S(t) +1) + R(t) 
=0 


from the SIR model (2) (see the left-hand sides of the SIR model equations and their 
sum). This implies that N(¢) can be assumed to be constant for the SIR model (2), i.e., 
N(t) = N(O) := N. With this assumption, we consider the state S(t) as N — I(t) — R(£). 
Then, we can describe the SIR model (2) as a two-dimensional model. 

Considering the magnitude of the numbers in Table 2, it is assumed that 


NZS. (4) 


Compared to the total size of the population, the number of outbreak cases is small. The 
N of South Korea is now known to be over 51 million. If the number of outbreak cases 
is much smaller than the total size of the population, the number N does not need to be 
exact to estimate system parameters of the SIR model (2) [21, 22]. 


From [16], 
B= KT 
= 


where T is the transmissibility of the infectious disease and x is the transmittable contact 
number by one infected subject per one unit time. 
Then, the model where the dimension is reduced by assumption (4) is 


T=xtl—vI, (5) 
R=vi. 

The initial state for the SIR model is provided by Table 2, given by 
[1(0), R(O)] = [2, 0]. (6) 


Now, we only need to know the parameter pair, (kT, v), to solve model (5) numerically. 
To solve the mathematical model and to evaluate the basic reproduction number of the 
model, we do not need to know each x and Tt necessarily, but we can use the value of the 
product, ie., KT. 

We search for the parameter pair (xt, v) such that can respond appropriately with the 
data in Table 2. To evaluate how closely system response is fitted to the data in Table 2, 
we employ a quantitative measure, the sum of the squared errors. Once we obtain the 
optimal values for the parameters with respect to this measure, we can estimate the 
basic reproduction number as described in the “Methods” section. 
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We define the measure /ff(-,-) as 


n 


felkiti wi) = .| > (Gi) — Ir)? + (Ri) - Rr(&)*) (7) 


k=0 


where m = 23 and k corresponds to the date. For example, k = 0 and k = 23 indicate the 
dates 20 May 2015 and 12 June 2015, respectively. J;(-) and R;(-) are from the simulation 
result of model (5) with the initial condition (6) and the corresponding parameters, i.e., 
the function arguments («;T;, vj). Rr(k) is the sum of the values of the “Deceased” and 
“Recovered” at k in Table 2. I7(k) is the value, Rr(k) subtracted from the value of the 
“Infected” at k. See Table 3. Based on the data of Table 2, I7(k) and Rr(k) along with k 
are shown in Table 3. 

To compare the quantitative measures for each pair of parameters, we consider the 
plane, i.e., 2-dimensional space, of the parameters, (xT, v). We can obtain a surface in 
3-dimensional space by plotting the corresponding measure as the value along the third 
axis. 

We explore the plane for wide ranges of parameters, and one of the results is shown in 
Fig. 1. This figure can help show the relation between f¢ and the system parameters. To 
effectively capture the characteristics of the measure, the ranges of «t and v in the figure 
are [0.00, 0.26] and [0.00, 0.12], respectively. 


Table 3 I;(k) and R7(k) of function (7), derived from Table 2 


k Ir (k) Rr(k) 
0 2 0 
] 3 0) 
2 3 0 
3 3 0) 
4 3 0) 
5 3 0 
6 5 0 
7 5 0 
8 7 0 
9 13 0) 
0 15 0 
1 18 0) 
2 24 ] 
3 29 1 
4 27 3 
5 32 4 
6 36 6 
7 58 6 
8 81 6 
9 86 9 
20 98 10 
21 109 13 
22 109 17 
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y oe KT 


Fig. 1 Plot of fe(«t, v) of (7) on the plane of (kt, v). fe(-, -) is a quantitative measure, the sum of the squared 
errors, that describes how the system response of model (5) with corresponding parameters is close to the 
data in Table 2. The ranges of «t and vin this figure are [0.00, 0.26] and [0.00, 0.12], respectively 


Based on characteristics shown in Fig. 1, we search for precise values for xt and v to 
minimize the measure, fe. Consider the parameter plane with the range [0.2050, 0.2150] 
of «t and with the range [0.0200, 0.0300] of v in Fig. 2. Then we can identify the pre- 
cise values for kt and v using Fig. 2. The system parameters are xt = 0.21153 and 
v = 0.026122 evaluated by Fig. 2 to minimise the quantitative measure f¢ 

The function f¢(-,-) describes the sum of squared error between the outbreak data 
of Table 2 and the simulation result of model (5) with the initial condition (6) and the 


0.215 


V 0.02 0.205 


KT 

Fig. 2 Plot of fe(«T, v) of (7) on the plane of (kt, v). Compared to Fig. 1, the ranges of «t and vin this figure 
are [0.2050, 0.2150] and [0.0200, 0.0300], respectively, and they can be used to obtain precise values of system 
parameters to minimize the quantitative measure fr 
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function arguments. Thus, if we find the parameters minimizing the function (7), then 
these parameters can be considered to correspond to the case of Table 2. 

Figure 3 shows the data presented in Table 2 and the state trajectories of model (5) 
with the parameter values obtained from Fig. 2. The top panel of Fig. 3 shows infected 
patients, and the bottom panel shows deceased or recovered patients. In both panels, the 
circle marks and the solid line display patient numbers based on Table 2 and the state 
transition of model (5), respectively. 

We derive the basic reproductive number for the SIR model (2) as 


Ro = — (8) 


in the “Methods” section, based on [16]. Equation (8) can be obtained alternatively by 
using the next generator approach [6, 8, 11, 27] to model (5). 

We determine the Ro of the MERS-CoV in South Korea in 2015 as 8.0977 (i.e., 
0.21153/0.026122). It is worth comparing with Ro of the 2014 Ebola outbreak, which had 
values of 1.5—2.5 [21]. 
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Fig. 3 Data presented in Table 2 and the state trajectories of model (5) with the estimated parameters. In 
the top panel, the circle marks plot the number of “Infected” patients from Table 2, and the solid line depicts 
the transition of state / of model (5). In the bottom panel, the circle marks indicate the sum of the numbers of 
“Deceased” and “Recovered” patients, and the solid line depicts the transition of state R of the model 
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Conclusions 

In this paper, we evaluated the basic reproduction number of the MERS-CoV outbreak 
that occurred in 2015 in South Korea, using officially reported data. We employed a 
mathematical dynamic model, the SIR model. We first fit the response of the SIR model 
to the epidemic curve data reported from the MERS outbreak. Then, we identified the 
system parameters of the model to estimate the basic reproduction number. 

Because there was no intervention to control the infection in the early phase of the 
outbreak, the data used here provide the best conditions to evaluate the epidemic char- 
acteristics of MERS, such as the basic reproduction number. An evaluation of Ro using 
epidemic data could be problematic if there are stochastic fluctuations in the early 
phase of the outbreak, or if the report is not accurate and there is bias in the data [11]. 
Such problems are not relevant to this study because the data used here were precisely 
reported and verified by [4]. 

We conclude this paper with the following discussion on future work to overcome the 


limitations of research, derived from assumptions in the paper. 


Further research direction 

¢ Behind the SIR model (2), there are several strong assumptions, one of which is a zero 
latent period, ie., the incubation period is zero. This implies that a patient becomes 
infectious immediately after infection. However, the incubation phase occurs dur- 
ing the course of the MERS outbreak. To address this weak point, in future work we 
could consider the 4-dimensional SEIR (i.e., susceptible-exposed-infectious-removed) 
model, which has been employed in [28] to study Ebola epidemic model. The addi- 
tional state in the SEIR model can help us deal with the latent period. 

¢ In this paper, we considered the epidemic curve data in [4] only from the early stage of 
the 2015 MERS outbreak in South Korea, where there was no intervention to control 
the spread. Accordingly, we evaluated Ro based on the data. In future work, we will 
also consider the epidemic data in [4] from the later (or closing) stage of the MERS 
spread in South Korea in 2015, so we can estimate the effective production number 
(ie., Reg), which is the production number resulting from interventions, such as edu- 
cation, quarantine, and the tracing of contacts by infected patients. By doing so, we 
can evaluate the effectiveness of each control measure on the spread of the infectious 
disease[29]. Eventually, such evaluation could help us improve public health policy. 
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